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We study bulk recombination in a medium with energetic disorder under continuous irradiation 
with excitation hght by using the multiple trapping (MT) model. Charge densities are shown to 
obey power law as a function of the intensity of excitation light, and the exponent is related to 
the dispersion parameter a of the MT model. The theory is applied to analyze charge densities in 
bulk heteroj unction photovoltaic cells measured by light-induced electron spin resonance (LESR). The 
', values of the exponent are consistent with those obtained from the transient decay of charge densities 
in time-resolved experiments. We show that the value of the exponent is robust against addition of 
^ . Gaussian trapping energy distribution which is due to shallow trap states. In order to corroborate 
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our theory, we also calculate the dependence of the open-circuit voltage on the light intensity on the 



CN ' basis of the same model and compare it with the corresponding experimental results. 
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I. INTRODUCTION 



In organic solar cells, it is important to suppress recombination between electrons and holes. 
So electron-hole recombination has been studied extensively in steady-state conditions in which 
the sample is continuously irradiated with an excitation light and in time-resolved conditions 
in which it is irradiated with a short pulse of the excitation light. 

If recombination between electrons and holes occurs by an ordinary bimolecular reaction, the 
carrier density in the case when the sample is continuously irradiated with the excitation light 
should be proportional to I^'^, where / is the light intensity. Recently Kuroda et al. (jl-3|) 
measured the dependence of the carrier density on the light intensity by using light-induced 
ESR (LESR). Although the observed data still have some scatter, they approximately followed 
a power law H, but the power 7 was much less than 0.5. Kuroda et al. explained this result 
by assuming that recombination occurred by a quadrimolecular reaction. 

By the same token, if recombination between electrons and holes occurs by an ordinary 
bimolecular reaction, the decay of carriers after they are generated by a short pulse of the 
excitation light should be proportional to where t is the time after the pulse. Recently 
several groups (j4-7|) measured the decay kinetics of carriers. The observed data followed a 
power law t"^ but (3 was much less than unity. This result was also explained by assuming that 
recombination occurred by multimolecular reaction with total reaction order larger than 2. 

In the present paper, we explain the above two types of experimental results on the basis 
of the multiple trapping (MT) model. The MT model has been used in order to elucidate 
the kinetics of charge recombination processes. js iSl The model used in the present paper 
was originally developed in 2002 in ref. 10|. In this paper the MT model which incorporates 
charge trap states with an exponential distribution of trapping energies was presented to an- 
alyze time-resolved kinetics of charge recombination in dye-sensitized solar cells and solved it 
analytically. In 2003 a similar model was used by Nelson to analyze time-resolved kinetics of 
charge recombination in organic solar cells, jol However, in Nelson's paper the model was not 
solved analytically but analyzed only by using numerical simulations. In 2009 we extended our 
model of 2002 to apply to time-resolved kinetics of charge recombination in organic solar cells, 
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and solved it analytically. 12|] We showed theoretically on the basis of this theory that the total 
density of charges decays with time as with /3 less than unity. In the MT model, charges 
perform transitions among trap states of different energies and the decelerated decay reflects 
the decrease in the detrapping rate as the distribution of charges am ong trap states of different 



s inside the band gap. 10Nl2| A similar explanation 



energies shifts toward deeper trap level 
was also given by other approaches. HQ 

If the total density of holes is changed in steady-state experiment by changing the light 
intensity, it changes distribution of holes among trap states with different trapping energies. 
This affects the open-circuit voltage (Vqc)- In order to corroborate our theory, in the present 
paper we also calculate the dependence of the open-circuit voltage on the light intensity on the 
basis of our model and compare it with the corresponding experimental results. 

In Sec. ini the MT model is generalized to apply to bulk recombination in the steady state 
conditions. Vqc is also calculated on the basis of the same model. Experimental results on 
steady-state and time-resolved conditions are analyzed in Sec. IIIII In Sec. IIVI the results are 
compared to those obtained by other models. The last section is devoted to the conclusion. 

II. THEORY 

A. Charge recombination dynamics 

Under the continuous irradiation with light, the excitons are generated. The electron-hole 
pairs are strongly bounded in excitons. The excitons can be dissociated into a hole and an 
electron at the interface between a donor phase and an acceptor phase in a bulk heteroj unction. 
The charge generation rate, G, is proportional to the light intensity, /, 

G = k,I. (1) 

The charges are initially produced in the free state. The motion of charges in bulk hetero- 
junction solar cells has been described by the MT model. JigI . IitI In the model, charges are 
trapped in the trap sites. We assume that if the detrapping rates of electrons and holes are 
different (This case is referred to as asymmetric case), carriers with the lower detrapping rate 



are always trapped but those with the higher detrapping rate are partly trapped and partly 
free. For the sake of convenience we hereafter assume in the paper that holes have the higher 
detrapping rate. The results for the opposite case are easily obtained from those presented. 
The case when the detrapping rates of electrons and holes are identical (This case is referred 
to as symmetric case) is analyzed in Appendix B. 

We assume that recombination between holes and electrons occurs once holes are detrapped. 
Therefore the recombination rate depends on the detrapping rates of holes. In time-resolved 
experiments distribution of holes among trap states with different trapping energies changes 



with time. This changes the detrapping rates of holes. We have shown theoretically in ref. [12 1 
that the total density of holes decays with time as with /3 less than unity instead of t~^. 
In this paper we analyze steady-state experiments. If the total density of holes is changed in 
steady-state experiment by changing the light intensity, it changes distribution of holes among 
trap states with different trapping energies. This changes the detrapping rates of holes. We 
show theoretically that the total density of holes depends on the light intensity as /''' with 7 
much less than 0.5 instead of I^'^. 

The trapping energy E of a. hole is measured relative to the HOMO of the donor. We 
denote the trapping energy distribution by g{E). We mainly treat the case when the trapping 
energy distribution is exponential. We also consider a more general case which includes an 
additional contribution of Gaussian trapping energy distribution. The detrapping rate constant 
is expressed as, 

kd{E) = iydexp{-E/kBT), (2) 

where is the detrapping frequency. We use to normalize rates and introduce the dimen- 
sionless energy defined by e = E/ksT. By thermal excitation, holes can be detrapped to the 
free state. The holes in the free state can either be trapped by vacant trap sites or recombine 
with trapped electrons. 

Let /(e) denotes the number density of holes trapped at trap sites with trapping energy e. 
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In the steady state the distribution /(e) obeys, 



kt{N -n) + krn •' ^ ^ 



oo 



A^^(e) -/(e) , _ 



' kt{N — n) + krU 

where the number density of trap sites is denoted by A^. n is the total number density of holes 
and related to the number density of holes trapped with energy e by, 

roo 

n= / f{e)de. (4) 
Jo 

The number density of electrons is the same as that of holes, kt and kj. are the rate constants of 
a hole in the free state for trapping and for recombination with a trapped electron, respectively. 
G, kt and fc^ are normalized by as mentioned above. The first term on the right-hand side 
represents the generation rate of holes with trapping energy e. The probability that a hole 
in the free state will be trapped in a trap state with trapping energy e is proportional to kt 
multiplied by the number density of the vacant trap sites given by Ng{e) — /(e), where Ng{e) 
is the total number density of trap states with trapping energy e. The probability that a hole 
in the free state will recombine with a trapped electron is proportional to the recombination 
rate multiplied by the number density of electrons, n. Therefore, the probability that a hole in 
the free state will be trapped into a site with trapping energy e can be written as, 

kt[Ng{e)-f{e)] ^ fc,[iV^(e) - /(e)] 

kt f^[Ng{e) - f{e)]de + Kn hN + {K - kt)n ^ ' 

This factor is multiplied by a hole generation rate in the free state, G, to obtain the generation 
rate of holes with trapping energy e. The second term on r.h.s. of Eq. (j3j) represents the 
decrease of holes with trapping energy e by detrapping. The physical meaning of the third 



term on r.h.s. was previously explained in detail. 

By integrating Eq. (|3]) over e, the fraction of trap sites occupied by holes, p = n/N, is shown 
to satisfy, 

POO 

l-p = y cieexp(-e)/(e)/iV, (6) 
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where g{e)de = 1 is used, y is given by 

If the carrier generation rate G is high, y is smalL In the steady state, the number density of 
holes trapped with trapping energy e is obtained by introducing Eq. ([6]) in Eq. ([3]) as, 

' l + '!/exp(-£) 

Equation ([8]) indicates that the number density of holes obeys a Fermi distribution as in equi- 
librium, and the quasi- Fermi energy of a hole relative to the HOMO of the donor can be defined 
by, 

eF = ln(|/). (9) 

Since y is given by Eq. (jTj), ei? depends on the fraction of trap sites occupied by holes, p = n/N , 
and the light intensity through G. By substituting Eq. ([7]) into Eq. ([9]), we find that the quasi- 
Fermi energy depends logarithmically on p/G, 

eF = Hkr/kt)+\n{Np/G). (10) 

In the following, we derive the relation between p and G. The density of occupied sites is 
given by Eq. (E]) regardless of the form of the trapping energy distribution. By substituting 
Eq. (IH]) in Eq. (jll) and noticing p = n/N, we obtain, 

p^fde , (11) 

Jo H-!/exp(-6) 

By expanding the integrand in terms of y, Eq. (ITT]) can be rewritten as, 

oo 

P = E(-^)"^(^)' (12) 

n=0 

where ^(n) denotes the Laplace transform of (7(e), ^(ra) = deexp{—ne)g{e). We note ^'(0) = 1 
from the normalization condition. 
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1. Case of low carrier generation rate 



a. Exponential trapping energy distribution Here we assume that g{E) decreases expo- 
nentially with increasing trapping energy, 

g{E) = {l/Eo)exp{-E/Eo), (13) 

where Eq is a parameter indicating characteristic trapping energy. The trapping energy distri- 
bution can be expressed in the dimensionless unit as, 

g{e) = aexp(— ae). (14) 

a = ksT/Eo is a key quantity in the MT model, a is called the dispersion parameter and 
given by the ratio of thermal energy to the characteristic trapping energy of trap sites. Small a 
compared with unity indicates that thermal energy is small compared with the typical trapping 
energy of trap states. Thermally assisted detrapping is suppressed when a is small. 

We obtain g{n) = a/{n + a) for the exponential trapping energy distribution. By using this 
equation, Eq, (fT2|) can be expressed as, 

p = 2^1 (l,tt, 1 + a, -y) , (15) 



where 2F1 (1, a, 1 + «, —y) is the hypergeometric function. [18| 

If y ^ 1, Eq. f|T5|) is simplified to the following equation by utilizing the asymptotic 
expansion of 2-^1 (1, a,l + a, —y), 



sm vra 



Equation ( IT6|) can be rewritten as a relation between the number density, p, and the charge 
generation rate, G, 



Since the charge generation rate is proportional to the light intensity (Eq. ([1])), the number 
density of holes obeys power law as a function of the light intensity and the exponent is given 
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by a/ (a + 1). Equation flTTI) is valid under the condition y ^ 1. By using Eq. ([7]) we can show 
that this condition is equivalent to 

7ra kr N , , 

sin(7raj /ct G 

The power law dependence of p on the excitation light intensity, Eq. (fT7|) . is related to the 
transient decay kinetics of the total number density of holes. When the light irradiation is 
switched off, the hole density decays to its value under the dark condition. We have analyzed 
such a transient kinetics. For the exponential trapping energy distribution, the transient time 



dependence of the hole density is derived as, 12 1 



Pit) - lA", (19) 

where a is again the dispersion parameter. The value of a in Eq. (fT9!) should be equal to that 
in Eq. (fTTl). 

b. Influence of shallow trap states We now consider a more general case which includes 
an additional contribution of Gaussian distribution in the trapping energy distribution, 

^(e) =paexp(-ae) + {1 - p) ^=1=— exp (- ) , (20) 

where p denotes the fraction of sites represented by the exponential trapping energy distribution. 
Cat is introduced to normalize the truncated Gaussian distribution in which e takes only positive 
values. It is given by, = (27rcr^) deexp [— (e — emp/ (2o"^)]. The density of states given 

by the sum of Gaussian and exponential distributions were successfully used to describe hole 
densities in blend films of RR-P3HT and PCBM. (loj If we substitute Eq. (jH]) together with 



Eq. (EOl) instead of Eq. in Eq. (ED, Eq. (US) is generahzed to. 



1 , ,^ exp 

p = p- ^+ ? — V 21 

sm vra y + exp {-emj 

7T<y 1 

~ P- + (1 - P) exp (e„) /y. (22) 

sm 7ra y°' 



In the derivation we used the following approximation, 

1 

t/-iexp(e) i/ + exp(e™) 



deg{e)—-^^ M ^ I y (23) 

1 + y ^exp e y + exp 
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and extended the lower limit of integration to — oo by assuming <^ a <^ and retained just 
the leading term by expanding l/[l + y~^exp(e)] about e = e^. We further took the asymptotic 
limit, y ^ exp (e^). By substituting Eq. ([7]) into Eq. (!22|) we obtain, 

P = P-. [tTt-] + 1-P exp e„ 24 

sm na \kr i\ p / kr i\ p 

When p(-7ra)/(sin7ra) > (1 — p) exp (e^) and/or in the limit of weak light intensity, Eq. ( p4|) 
leads to the same power law as that given by Eq. f lT7|) . As long as the tail of trap states can 
be approximated by an exponential trapping energy distribution (Eq. (fT4|) ). the exponent for 
the power law dependence of the hole density on the continuous light intensity is robust and is 
given by a/(l + a) even in the presence of additional shallow states. 

In the ordinary theory of bimolecular recombination in the absence of trap states, the hole 
density as a function of continuous light intensity is given by power law with the exponent 0.5. 
The same result is obtained when we consider Gaussian trapping energy distribution, 

with a small width of energy distribution, a ^ e^- By setting p = in Eq. (!24|) we get 



ktG 



p = exp(e^/2)y^-, (26) 
as expected. Note that Eq. fl26|) is obtained by assuming y ^ exp (e^) ^ 1- 



2. Case of high carrier generation rate 

If y ^ 1, we obtain the following equation by retaining the first two terms in Eq. (lT2l) and 
further using Eq. (j7]), 

p ^ 1 + . (27) 



The fraction p of trap sites occupied by carriers saturates when the light intensity is high. 
The saturation results from the filling of the trap levels inside the band gap up to the shallow 
trap levels. Note that Eq. ( 1271) is derived without assuming any form of the trapping energy 
distribution. 
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B. Open-circuit voltage 



According to the MT model, distribution of holes among trap states with different trapping 
energies can be characterized by the quasi- Fermi energy. Figure [1] shows the hole density as 
a function of the trapping energy for the exponential (A) and Gaussian (B) trapping energy 
distributions. In both (A) and (B), the hole density is shown for two light intensities which are 
different by a factor of 10. Suppose that the quasi-Fermi energy at the low light intensity is 
the same for both (A) and (B). The thick dashed line in (A) and (B) indicates the quasi-Fermi 
energy when the light intensity is increased 10 times. For this change of the light intensity the 
quasi-Fermi energy for the exponential trapping energy distribution decreases more than that 
for the Gaussian trapping energy distribution. In other words, the amount of change in the 
quasi-Fermi energy by increasing light intensity depends on the form of g{e). It implies that the 
open-circuit voltage estimated from the quasi-Fermi energy is affected by the trapping energy 
distribution and the light intensity; the open-circuit voltage in the case of the exponential 
trapping energy distribution depends on the light intensity more strongly compared to that in 
the case of the Gaussian trapping energy distribution. 

We now study quantitatively the influence of the trapping energy distribution on the open- 
circuit voltage under the illumination of continuous light. In open-circuit condition, the current 
is not extracted and the photogenerated carriers in principle disappear only by recombination. 
When recombination is in equilibrium with generation of charge pairs, the open-circuit voltage 
is given by the difference between the quasi-Fermi energy of a hole and that of an electron both 



measured from the same level, for example, the HOMO level of the donor, 20|, [21 1 



eVoc = [ksT) (e^ - ep^ , (2^ 



where e^^ is the quasi-Fermi level of an electron measured from the HOMO level of the donor. 
We first consider an asymmetric case and assume that the detrapping frequency of electrons, 
is smaller than that of holes, u^, ^ = vf^ jvd ^ 1- When ^ = v"^^ jvd ^ Ij the quasi-Fermi 
energy of a hole can be approximated by Eq. ffTOj) . In Appendix A, the quasi-Fermi energy 
of an electron in this case is calculated as Eq. (]A.6p . This quasi-Fermi energy of an electron, 
e'^'', is measured from the LUMO level of the acceptor in the downward direction. (See Fig. [2] 
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FIG. 1: (A) Hole density as a function of the trapping energy for the exponential trapping energy 
distribution with a = 0.2. Dotted line represents the exponential trapping energy distribution of trap 
states given by Eq. (jl4p . The thin solid line represents the distribution given by f{e)/N with the 
quasi- Fermi energy equal to 0.173 [eV] which is indicated by the thin dashed line. The thick solid 
line shows the distribution when the light intensity is increased 10 times. The quasi-Fermi energy 
for this case is equal to 0.125 [eV] and is indicated by the thick dashed line. (B) Hole density as a 
function of the trapping energy for the Gaussian trapping energy distribution with the mean and the 
variance given by 0.1 [eV] and 0.03 [eV], respectively. Dotted line represents the Gaussian trapping 
energy distribution of trap states given by Eq. (j25p . The thin solid line represents the distribution 
given by f{e)/N with the quasi-Fermi energy equal to 0.173 [eV] which is indicated by the thin dashed 
line. The thick solid line shows the distribution when the light intensity is increased 10 times. The 
quasi-Fermi energy for this case is equal to 0.144 [eV] and is indicated by the thick dashed line. 



for the energy coordinates.) The quasi-Fermi energy of an electron measured from the HOMO 
level of the donor in the upward direction is given by, 



P - acceptor LUMO level - donor HOMO level - e"-p' 

acceptor LUMO level — donor HOMO level — log {^kf/kr) — ep, 



(29) 
(30) 



where e'^^ given by Eq. f lA.6P is substituted. We treat the two cases when the trapping 
energy distribution of holes is Gaussian and when the trapping energy distribution of holes is 
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FIG. 2: Schematic representation of e and e' coordinates of trap states. 



exponential. 



1. Case of low carrier generation rate 

a. Gaussian trapping energy distribution When the trapping energy distribution is Gaus- 
sian, the quasi- Fermi energy of a hole measured from the HOMO level of the donor is obtained 
by substituting Eq. ([26D in Eq. ([TOD, 

= (1/2) [e„ + \n{kr/kt) - ln(G/iV)] . (31) 



By substituting Eqs. ( l30l) and ( ISTl) in Eq. ( !28l) . Vqc is obtained as a linear function of the 
logarithm of light intensity, /, 



Voc = ci + {kBT/e) In/, 



(32) 



where Eq. ([T]) was used and ci is a constant independent of the light intensity. According to 
Eq. ( l32l) . the relation between Vqc and In J is linear with a positive slope equal to ksT /e. This 



is indeed observed experimentally. 



22 



23| Equation fl32|) was also derived by assuming a single 



trap state both for a hole and an electron. 22h24| 
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b. Exponential trapping energy distribution When the trapping energy distribution is ex- 
ponential, the quasi- Fermi energy of a hole is obtained by substituting Eq. (fT71) in Eq. (fTOl) . 



a 



In 



ira 



sm na 



+ \n{K/kt) - \n{G/N) 



By substituting Eqs. fl5U|) and ( 1551) in Eq. (1251) Vqc is obtained as, 

Voc = C2 + [2/(« + 1)] {kBT/e) In J, 



(33) 



(34) 



where C2 is a constant independent of the light intensity. The slope of Vqc against In / could be 
in the range between kBT/e and 2kBT/e in this case. As already mentioned, Eq. (1T7|) is valid 
for [na/ sinln a)]{kr/kt){N/G) ^ 1. Therefore, Eq. is also valid under this condition. The 
slope larger than kBT/e is recently reported in bulk hetero junction solar cells, where charge 



transport is trap-limited. 



22 
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26| 



So far, we have assumed that the hopping frequency of holes is larger than that of electrons. 
In the opposite situation, a in Eq. (134|) should be replaced by that for electrons. 

2. Case of high carrier generation rate 

By using Eq. ([7]) we can show that the condition of j/ <^ 1 is equivalent to 

na kr N 



sin(7ra) kf G 

By using Eq. (l27jl . Eq. ((7]) can be expressed as. 



< 1. 



(35) 



y 



m 

hG' 



hG_ 

krN 



(36) 



where Eq. ( l35l) and ^'(1) < 1 are used. As a result, = Iny is obtained as, 

eF^ln{kr/kt)-ln{G/N). (37) 

By substituting Eqs. ( 130|) and ( 1371) in Eq. ( l28l) . Vqc is obtained as a linear function of the 
logarithm of light intensity, /, 



Voc = ci + 2{kBT/e)\nI. 



{31 
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It should be reminded that when the hght intensity is high, this relation is derived without 
assuming any specific forms of 5'(e). The slope of Vqc vs ln(/) is given by 2kBT/e regardless 
of the forms of g{e) as a result of trap filling effect explained below. 



3. Case of identical kinetic parameters of electrons and holes 

Recently, the relations given by Eq. fl32|) and Eq. flMj) have also been found by numerical 
simulations when the kinetic parameters of electrons and holes are identical. 27|] Equation 
(132|) is found for the Gaussian trapping energy distribution and Eq. fl3ll) is found for the 
exponential trapping energy distribution, j^^l In the simulation, the drift-diffusion equation 
was solved together with the Poisson's equation, and the carrier density in the trap states, pt, 
was approximated as pt ~ exp (^—E'j/Eoj in terms of the quasi-electrochemical potential E'jr 
determined self-consistently. The results were interpreted by using pt given above with spatially 



homogeneous |27l. |28| 



Here, we calculate Vqc analytically in the case when the kinetic parameters of electrons and 
holes are identical. If y given by Eq. ( ]B.4[) is larger than unity, the quasi-Fermi energy is given 
by Eq. (IB. 81) . Since we have e'p^ = ^f, Vqc is obtained as Eq. (IMl) by using Eq. (IB. 81) and 
Eqs. ( I28l) -(l29l). Equation ( 134|) is valid when y given by Eq. ( ]B.4[) is larger than unity which is 
equivalent to 

2ktkr N , , 

ftWG » ' (^^) 

Let us consider the opposite limit y ^ 1 which is equivalent to, 

2ktkr N , , 

This condition is satisfied under the strong light intensity leading to high carrier generation 
rate, G. The quasi-Fermi energy is given by Eq. (IB.lip . Since we have = ep, Vqc is 
obtained as Eq. §^ by using Eq. (IRTTI) and Eqs. (|28D-(|29D. 
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III. ANALYSIS OF EXPERIMENTAL RESULTS 



A. Dependence of stationary hole density on excitation light intensity 

Charge carrier density in the composite of regioregular poly3-hexyhhiophene (RR-P3HT) 
and 6,6-phenyl-C61-butyric acid methyl ester (PCBM) was recently measured by LESR under 
the continuous irradiation with light, jlj] LESR technique has been applied to probe the holes 
on the polymer and fullerene radical anions. jl-3| 

In Fig. [3l the hole density measured by LESR is shown against excitation light intensity for 
different PCBM concentrations. [l| We analyzed these results by using power law dependence 
on the excitation light intensity. As shown in Fig. [3l the values of a estimated by using Eq. 
(frfl) are within the range of 0.26 — 0.35. These values can be assessed by analyzing the transient 
kinetics. 



B. Time-resolved decay of hole density 

As shown in Fig. HI the decay of the hole density in time-resolved experiments obeys power 
law and we obtained a as 0.3 for 5 mol % of PCBM by using Eq. ( |T9l) . This value is consistent 
with a = 0.31 estimated from Fig. [3] by using Eq. f|T71) . The difference between the fit with 
a = 0.31 and the best fit with a = 0.3 is small as shown in Fig. HI The results indicate the 
applicability of MT model to interpret the experimental data on both the steady state and time- 
resolved conditions. Unfortunately, the experimental results shown in Fig. [3] are not sufficiently 
accurate. More accurate and systematic measurements of experimental data are badly needed 
in order to make a more quantitative comparison between our theory and experiment. 

Now we analyze the temperature dependence of the decay kinetics. Unfortunately, detailed 
data on the temperature dependence of the decay kinetics are not available for this sample. So, 
we analyze the temperature dependence of the decay kinetics for other samples. 

Figure [5] (A) shows the decay of the hole density for RR-poly(3- 
octylthiophene) (P30T) /fullerene at temperatures, 60, 80, and 100 K. Q All results 
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FIG. 3: Hole density against excitation light intensity of 670 nm. RR-P3HT/PCBM for different 
PCBM concentrations measured at 60 K. [l| The lines indicate the fitting by Eq. (jl7p . a values 
extracted from the fitting are a = 0.35 for 55 mol %, a = 0.26 for 36 mol % and a = 0.31 for 5 mol 
%. For 5 mol % the value of the exponent is a/(a + 1) = 0.236. 

indicate power law decay at long times. The value of the exponent a obtained from the decay 
kinetics is shown in Fig. |5] (B) as a function of temperature, a can be well approximated by a 
linear function of temperature and is almost proportional to temperature. This is consistent 
with the MT model which incorporates an exponential trapping energy distribution. It predicts 
a = ksT / Eq. Power law decay of the hole density and linear temperature dependence of a 
is also found for RR-P3HT/fullerene as shown in Fig. [61 The dispersion parameter a should 
be proportional to the temperature if the trapping energy distribution is strictly exponential. 
The deviation indicates that the trapping energy distribution is not strictly exponential. 
In Figs. [5] and |6] the same acceptor was used, but the length of alkyl chains was changed. 
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FIG. 4: Normalized hole density against time after the switch-off of the excitation light of 670 nm 
with 2.7 mW/cm2. RR-P3HT/PCBM with 5 mol % PCBM measured at 60 K. Q The sohd line 
indicates the best fit by power law kinetics with a = 0.3. The dashed line shows the fitting with 
a = 0.31 obtained by analyzing Fig. [3l 

If experimental data are more accurate and systematic, more information will be obtained on 
the effect of alkyl chain length on the trapping energy distribution by comparing Figs. [5] and 
El More accurate and systematic measurements of experimental data are badly needed in order 
to make a more quantitative comparison between our theory and experiment. 
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FIG. 5: (A) Normalized hole density against time. RR-P30T/fullerene with 5 mol % fullerene 
measured at 60, 80, and 100 K. ^ The excitation light of 700 nm with 1 mW/cm^ is switched off at 
time zero. Squares, diamonds, and circles correspond to 60K, 80K, and lOOK, respectively. The solid 
lines are the fits by power law kinetics, Eq. (|19p . (B) a values determined by fitting to Eq. (|19|) as a 
function of temperature. The line represents the result of fitting by a = 0.01 + T/333. 
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FIG. 6: (A) Normalized hole density against time. RR-P3HT/fullerene with 5 mol % fullerene 
measured at 40, 60, and 80 K. [2] The excitation light of 700 nm with 0.25 mW/cm^ is switched off at 
time zero. Squares, diamonds, and circles correspond to 40K, 60K, and 80K, respectively. The solid 
lines are the fits by power law kinetics, Eq. (fT9]) . (B) a values determined by fitting to Eq. (fT9]) as a 
function of temperature. The line represents the result of fitting by a = 0.162 + T/444. 
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IV. DISCUSSION 



We have assumed that recombination between holes and electrons occurs once holes are 
detrapped. Therefore the recombination rate depends on the detrapping rates of holes. In time- 
resolved experiments distribution of holes among trap states with different trapping energies 
changes with time. This changes the detrapping rates of holes. We have shown theoretically 
that the total density of holes decays with time as with /3 less than unity instead of t~^. If 
the total density of holes is changed in steady-state experiment by changing the light intensity, it 
changes distribution of holes among trap states with different trapping energies. This changes 
the detrapping rates of holes. We have shown theoretically that the total density of holes 
depends on the light intensity as /''' with 7 much less than 0.5 instead of P'^. 

In the multimolecular reaction model the recombination between electrons and holes is as- 
sumed to occur by multimolecular reaction between electrons and holes. When total reaction 
order is m, the recombination rate is proportional to the mth power of the total density of holes, 
if the total densities of holes and electrons are equal. In the multimolecular reaction model 
the recombination rate is assumed to be proportional to the densities of holes and electrons 
and the rate coefficient is implicitly assumed to be a constant. If this assumption is made, 
in steady-state experiments the total density of holes depends on the light intensity as Ji/"*. 
However, it should be pointed out that the rate coefficient in this case is not a constant for the 
following reason. If the density of holes is changed in steady-state experiments by changing the 
light intensity, it changes the distribution of holes among trap states with different trapping 
energies. This changes detrapping rates of holes and therefore changes recombination rates 
with electrons. In our opinion this is the reason the reaction order larger than 2 is obtained, if 
experiments are analyzed on the basis of the multimolecular reaction model. The same criti- 
cism also applies to the analysis of time-resolved experiments on the basis of the multimolecular 
reaction model. 

Previously, the trap-limited recombination has been interpreted by using the Shockley-Read- 



Hall(SRH) theory. 29| The original SRH theory of recombination assumes a single trap state 



in a material composed of a single phase. The original SRH model can be extended to the case 



19 



28 



However, it should be 



of trap states with an exponential distribution of trapping energies, 
pointed out that since the original SRH theory assumes the steady state, it cannot be applied 
to analyze transient kinetics. 

There is another important assumption in the SRH model. It assumes that a trap state, say, 
for an electron, can trap an electron efficiently, if it is empty, not only from the conduction band 
but also from the valence band. In our opinion, this assumption is reasonable for inorganic 
semiconductors but not for organic semiconductors. In organic semiconductors the electron 
trapping energy distribution does not extend very deep inside the band gap from the conduction 
band edge. In other words, trap states are located close to the conduction band edge but far 
away from the valence band edge. As a result, a trap state can trap an electron efficiently from 
the conduction band but not efficiently from the valence band. In our model we assume that 
trapping and detrapping of a hole from/to the valence band occur efficiently and recombination 
occurs efficiently between free holes and trapped electrons. 



There are several differences between our theory and those in refs [27[ and [28|. Differences 
are as follows. 

In our model the kinetic parameters of holes and electrons are in general different, while in 
their models it is assumed that electrons and holes have identical kinetic parameters. 

In our model the spatial variation of electrostatic potential is not taken into account, while 
in their models it is taken into account. When the spatial variation of the carrier density is 
small, their results agree with ours. 

Our model is solved analytically, while their models are analyzed mainly by using numerical 
simulations. In our theory the carrier distribution among the trap states are calculated rigor- 
ously, while in their models it is calculated only approximately. In this paper, we have solved 
the symmetric case as well analytically to compare with their corresponding results obtained 
numerically and with our results on the asymmetric case. 

In their models, the direct recombination between free holes and free electrons was also 
taken into account. 



27 



It was shown that the open circuit voltage can be limited by direct 



recombination when the trapping rates of free carriers are very small. 



Direct recombination 



between trapped holes and trapped electrons was also taken into account in a simulation, though 
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recombination between free carriers and trapped carriers of opposite charge was considered to 
be dominant. In this paper we assume that recombination occurs between free holes and 
trapped electrons. We assume that holes are mostly in trap states and can recombine with 
trapped electrons once they are detrapped to the conduction band and become free. We do 
not think that in organic semiconductors holes trapped near the valence band edge recombine 
efficiently with electrons trapped near the conduction band edge. 

It should also be pointed out that when they derived Eq. for the exponential trapping 
energy distribution, they implicitly assumed that the light intensity is low, so that the genera- 
tion rate is also low. In this case the condition given by Eq. fl39|) holds and Vqc is obtained as 
Eq. (134|) . However, in the opposite limit given by Eq. (140|) we have shown that Vqc is obtained 
as Eq. f p8|) as already shown. 

The above results are explained on the basis of the trap filling effect. Trap states have 
different trapping energies. When the light intensity is low, only deep trap levels inside the 
band gap are occupied by charges. When the light intensity is high, shallow trap levels close 
to the conduction band are also occupied. Charges occupied in shallow trap states are easily 
detrapped to the conduction band, so they behave like quasifree charges. Accordingly the slope 
of Vqc vs In/ becomes independent of the form of g{e). The trap filling effect is important 
when the light intensity is high and the number of trap states is small. On the other hand, the 
condition given by Eq. (l39l) holds when the light intensity is low and the number of trap sites 
is large. In this paper we calculate Vqc in both limits of (a) low light intensity with large trap 
site number and (b) high light intensity with small trap site number. 

V. CONCLUSION 

In this paper, we have analyzed recently observed power law dependence of the charge 
density on the intensity of excitation light by using the MT model. The observed relation 
between the hole density n and the excitation light / obeys power law, n P, where 7 is 
smaller than 0.5. Previously, the exponent 7 smaller than 0.5 was explained by assuming 
quadrimolecular recombination among holes and electrons. By using the MT model, we have 
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successfully explained the experimental results of 7 < 0.5. To be more precise, the exponent is 
given in terms of the dispersion parameter a as 7 = a/(l + a) and a should be smaller than 1 
when the deep trap levels inside the band gap is described by an exponential trapping energy 
distribution. 

The value of the exponent is consistent with that obtained from the independent experi- 
ment on the transient decay of charge density. Analysis of the temperature dependence of the 
exponent obtained from the transient decay shows that the temperature dependence of the ex- 
ponent is consistent with the MT model. The value of the exponent is robust against addition 
of Gaussian trapping energy distribution to the usual exponential trapping energy distribution, 
the former representing shallow trap states and the latter representing deep trap states. Even 
under the presence of shallow trap states the exponent characterizes the exponential trapping 
energy distribution of deep trap states. 

We also show that Vbc depends linearly on the logarithm of light intensity and the slope 
is equal to /c^T/e when only the shallow trap states are present. Otherwise, the same linear 
relation holds but the slope is larger than /c^T/e. 
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Appendix A. The detrapping rates of electrons and holes are different (Asymmetric 
case) 



We consider the case when the detrapping frequency of electrons is smaller than that of 
holes, namely ^ = v^f^ jvd ^ 1- Let /e(e') denote the number density of electrons trapped at 
trap sites with the trapping energy e'. The trapping energy e' of an electron is measured relative 
to the LUMO of the acceptor. In the steady state, /e(e') obeys the equation similar to Eq. (|3]), 



kf\N,-n) + 0^n kt{N-n) + krn 

where the number density of trap sites for electrons is denoted by Nf,, the second term on the 
right-hand side represents loss of trapped electrons with trapping energy e' due to recombination 
with generated holes, and the last term represents loss of trapped electrons with trapping energy 
e' due to recombination with holes detrapped into the free state. 

When e e-'' fe{e')de' < e-'f{e)de, ki'^ < kl"\ and K < h, Eq. (jAl]) simplifies into, 

c k 1 r°° 

Ggeie') = ie~'Ue') + —Ue') + fe{e')^^j^ | e~^f{e)de. (A.2) 
If the fraction of trap sites occupied by holes is small, Eq. (E]) can be rewritten as. 

By substituting Eq. flA.SP into Eq. (IA.2p . we obtain, 

G 1 + i/eexp (-e') 



when p Ne/N, where jje is given by. 



Ve = -^Pi. (A.5) 



Uf. can be further rewritten as Ue = ^{kt/kr)y in terms of y defined by Eq. ([7]). Therefore, the 
quasi- Fermi energy of an electron relative to the LUMO of the acceptor is given by, 

e'J:^ = log = log {ih/K) + ep. (A.6) 
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Appendix B. The detrapping rates of electrons and holes are identical (Symmetric 
case) 



We first note that Eq. (lA.ip holds without any restriction on the kinetic parameters of 
electrons and holes. If the kinetic parameters of electrons and holes are identical, Eq. (lA.ip 
can be simplified as, 



^ ^ [Ng{e)-f{e)]-krf{e) _ ^ h [Ng{e) - /(e)] - Kf{e) 



kt{N — n) + krU Jq 



kt{N — n) + krU 
By integrating Eq. (IB.ip over e, we obtain, 

^_ 2k,.n (G + /-e-7(6)&) 
kt{N — n) + krU 

Equation (]B.1|) can be rewritten by using Eq. (]B.2p as 



e-7(e)rfe. (B.l) 



(B.2) 



Ng{e) 



kt + kr J 1 + ?/exp(-e) 



(B.3) 



where y is given by. 



y 



N 2kr 



-P- 



(B.4) 



G kt + kr 

It is interesting to note that if we introduce ep through y = exp^ep), Eq. ( IB. 31) represents the 
Fermi-Dirac distribution with the quasi- Fermi energy given hj ep. By introducing Eq. flB.3|) 
into Eq. (j4]), we obtain, 

kt 



P 



kt -\- kr Jq 



de- 



9[^) 
1 + yexp{-e) 



(B.5) 



Equation (IB.SP is expressed as. 



kt 



P 



-2-Fi (1,Q!, 1 + a, -y) 



kt ~\~ kj- 

when the trapping energy distribution is exponential, Eq. ( fT4l) . 

Following the derivation of Eq. f[T7|) from Eq. (IT^ we have from Eq. fIB.Sp for ?/ ^ 1 

Txa kt 



(B.6) 



P 



sinfvra) kt + K 



2fc^ N 



(B.7) 



24 



Equation flB.7p shows that the number density of holes obeys essentially the same power law 
dependence on the light intensity as shown in Eq. (1T71) . By introducing Eq. ( IB .71) into Eq. 
(lB.4p . the quasi- Fermi energy of a hole given by ei;' = ln(?/) is expressed as, 



a 



In 



sinvra 



In 



2 /l^ \z^. 



- ln(G/iV) 



(B.8) 



While following the derivation of Eq. (!27|) from Eq. Eq. (ITT]) we have from Eq. (IB.SP for 

y > 1, 



p 



kt + kr 2krg{l) N 



(B.9) 



kt + K G ^ 

In this limit, the number density saturates as the generation rate, G, increases. The fraction 
p of trap sites occupied by carriers saturates when the light intensity is high. The saturation 
results from the trap filling of the deep trap levels inside the band gap up to the shallow trap 
levels. Note that Eq. ( 1B.9|) is derived without assuming any form of the trapping energy 
distribution. ^(1) = a/{l + a) is obtained for the exponential trapping energy distribution. By 
using Eq. (IB. 91) . Eq. (]B.4p can be expressed as, 



^(1) + 
/by' 



{kt + Kf G 



n -1 



2 Jxirp 

N 



N 



{kt + KYG' 

where Eq. pUj) and 1 > ^(1) are used. As a result, ep = In?/ is obtained as, 

2 }\if 



(B.IO) 



ep — In 



{kt 



kr)"^ 



-\n{G/N). 
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